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Abstract. We describe a new kludge scheme to model the dynamics of generic extreme-mass- 
ratio inspirals (EMRIs; stellar compact objects spiraling into a spinning supermassive black hole) 
and their gravitational-wave emission. The Chimera scheme is a hybrid method that combines 
tools from different approximation techniques in General Relativity: (i) A multipolar, post- 
Minkowskian expansion for the far-zone metric perturbation (the gravitational waveforms) and 
for the local prescription of the self- force; (ii) a post-Newtonian expansion for the computation 
of the multipole moments in terms of the trajectories; and (iii) a BH perturbation theory 
expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesies. The 
EMRI trajectory is made out of Kerr geodesic fragments joined via the method of osculating 
elements as dictated by the multipolar post-Minkowskian radiation-reaction prescription. We 
implemented the proper coordinate mapping between Boyer-Lindquist coordinates, associated 
with the Kerr geodesies, and harmonic coordinates, associated with the multipolar post- 
Minkowskian decomposition. The Chimera scheme is thus a combination of approximations 
that can be used to model generic inspirals of systems with extreme to intermediate mass 
ratios, and hence, it can provide valuable information for future space-based gravitational-wave 
observatories, like LISA, and even for advanced ground detectors. The local character in time of 
our multipolar post-Minkowskian self-force makes this scheme amenable to study the possible 
appearance of transient resonances in generic inspirals. 



1. Introduction 

The modeling of gravitational wave (GW) sources is key for the success of GW astronomy. 
Detectors require precise theoretical waveforms to extract and characterize signals buried in the 
noise. One of the main GW sources for space-based observatories, like the Laser Interferometer 
Space Antenna (LISA) [H(2]) that requires accurate templates for detection and analysis are 
extreme-mass-ratio inspirals (EMRIs) [3j. These events consist of a small compact object (SCO), 
such as a stellar-mass black hole (BH) or a neutron star, spiraling in a generic orbit into a 
spinning, (super)massive black hole (MBH), in the regime where GW emission is relevant. In 
such inspirals, the SCO spends up to millions of cycles in close orbits around the MBH, possibly 
with large pericenter velocities and eccentricities, sampling the strong gravitational field of the 
MBH. 



The relevant EMRI phase for data analysis is the inspiral, during which the number of 
cycles accumulated scales with the inverse of the mass ratio q = m±/M, [0(1O _4 )-0(1O~ 6 )]. 
During EMRI plunge and merger, the accumulated number of cycles scales with the MBH mass 
(M.) only. Moreover, the signal-to- noise ratio (SNR) is in the range O(10)-O(10 2 ) for EMRIs 
at realistic distances, and it scales linearly with the total number of cycles. Therefore, the 
contribution of the merger/plunge phase to the SNR is a fraction 0(m*/M,) relative to the 
total inspiral contribution. In addition, the EMRI ringdown phase is not detectable, as the SCO 
barely perturbs the MBH geometry as it plunges and crosses the MBH's event horizon. 

As a consequence of the large number of GW cycles contained in EMRIs, they can provide a 
detailed map of the MBH geometry, allowing us to extract the EMRI physical parameters with 
high precision [3] . This information is crucial to test the Kerr hypothesis J5j |6l El HH] arid 
constrain modified gravity theories (see, e.g. [TTJ, [12l [THl [13] ) . EMRI observations will also have 
an important impact on astrophysics (for a review see [3]). 

The extreme mass ratios involved in the problem also lead to the appearance of two different 
spatial and temporal scales. The two different spatial scales are due to the different MBH and 
SCO sizes, m*/M. <C 1. The two different temporal scales are the orbital one(s) and the one 
associated with radiation reaction, 7o rbital /T RB ~ m*/M, <C 1. Due to this, a fully relativistic 
numerical evolution of EMRIs is currently unfeasible. 

A framework to model EMRIs that exploits the extreme mass ratios involved is BH 
perturbation theory, where one treats the SCO as a small perturbation of the MBH background 
geometry. In this context, the inspiral is modeled as due to a local self-force. This force 
is composed of the metric perturbations generated by the SCO, regularized by eliminating 
divergences due to the SCO's point particle description. The SCO's motion is then governed 
by the MiSaTaQuWa equation of motion, derived in |14| [T5] . This equation is the foundation 
of a self-consistent scheme to describe EMRIs in an adiabatic way, by coupling this equation 
of motion to the partial differential equations that describe the perturbations produced by the 
SCO. 

At present, the gravitational self- force has been computed for the case of a non-rotating MBH 
using time-domain techniques [T7] . Progress is currently being made toward the calculation 
of more astrophysically relevant (spinning) EMRIs |18| . Given the number of cycles contained 
produced by an EMRI in one year and the present complexity of self-force calculations, we cannot 
expect to generate complete waveform template banks by means of full self-force calculations. 
Instead, it has been proposed that the goal of these studies should be to understand all the details 
of the structure of the self- force, so that one can formulate efficient and precise algorithms 
(phenomenological models) to create the waveforms needed for LISA data analysis, perhaps 
complementing existent approximation schemes. 

In this paper, we summarize the results of [19] where a new approximate scheme was 
proposed to model EMRIs: the Chimera framework. This scheme can be thought of as a 
new kludge model that combines different approximation methods, including BH perturbation 
theory, the multipolar post-Minkowskian formalism of Blanchet and Damour [20] and post- 
Newtonian theory. The Chimera scheme could be used for qualitative descoping studies, i.e. to 
determine the accuracy to which parameters could be extracted given an eLISA EMRI detection 
as a function of detector properties. The Chimera scheme could also be used to study how 
certain resonances, found by Flanagan and Hinderer |21| . affect such parameter estimation 
and detection. Chimera waveforms, however, would require improvements before they can be 
used as realistic templates in data analysis. The Chimera scheme is indeed amenable to such 
improvements, which will be the focus of future work. 



2. Ingredients of the Chimera Scheme 

2.1. MBH Geometry and Geodesic Motion 

The geometry of the MBH is described by the Kerr metric. In Boyer-Lindquist (BL) coordinates, 
(^bl) = r ) 0, (j)), the line element is given by 
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with / = 1 — 2M,/r. The function 
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A has two roots: r_j_ = M, ± yM, 2 — a 2 ; the positive one r , 
of the event horizon. The Kerr geometry is a stationary, axisymmetric and vacuum solution to 
the Einstein equations, which has has an additional symmetry described by a Killing tensor. 

In the limit where the SCO has zero mass (m t — > 0), it follows timelike geodesies of the MBH 
geometry which, in terms of proper time, can be parameterized as z^{t) = (i(r), r(r), 9(t), 4>(t)). 
Given the symmetries of the Kerr metric, each of these geodesies can be fully characterized by 
three constants of motion (per SCO unit mass): The energy (stationarity) , E; the angular 
momentum in the MBH spin direction (axisymmetry), L z , and the Carter constant (Killing 
tensor symmetry), Q (or C = Q — (L z — aE) 2 ). Using these symmetries we can completely 
separate the geodesic equations so that they can be written in the form: 
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where S 2 = (r 2 + a 2 ) 2 - a 2 A sin 9 . For the purposes of EMRI modeling we restrict ourselves 
to bound geodesies (orbits), which can also be characterized in terms of orbital elements: the 
eccentricity e, the semilatus rectum p, and the inclination 9 inc . The first two are defined in 
terms of the turning points (extrema) of the radial motion, the pericenter and apocenter (r pcri 
and r apo ), by the Newtonian relations: r pori = pM,/ (1 + e) and r apo = pM 9 / (1 — e) . Likewise, the 
inclination angle 9 inc is associated with the turning point of the polar motion, i.e. the minimum 
of 9 (9 min € [0, 7r/2]), via the relation: 9 inc = sign(L 2 ) (vr/2 — min ) . Another useful definition of 
the orbital inclination angle is: cost = L Z /\J L\ + C . 

Alternatively, the orbits can also be characterized in terms of three fundamental frequencies 
(see e.g. [221 E3 El]) with respect to the BL time t (they can be also constructed using proper 
time or any other time): f2 r , associated with the radial motion (from periapsis to apoapsis 
and back); fig, associated with polar motion; and O^, associated with azimuthal motion. 
Precessional orbital effects are due to mismatches between these frequencies, and they can be 
used to decompose, among other things, the GWs in a Fourier expansion [23J. 



2.2. Chimera Osculating Trajectories 

The foundations of the gravitational backreaction of a moving massive body on its own trajectory 
were laid down in the papers by Mino, Tanaka, and Sasaki |14j and Quinn and Wald [15]. These 
works established that the trajectory of a point mass can be understood as geodesic motion in 
the background geometry modified by a local force, the self-force, given in terms of the gradients 



of the regularized^ first-order metric perturbations. Then, the equation of motion for the SCO, 
the MiSaTaQuWa equation, is 



d 2 z 



+ T^u^ = m- 1 F s a F , (6) 



dr 2 ^ 
where the self- force, F® p , is given by 

F s % = ~\m* (g aX + u a u x ) u*V (2V„/& - V x h%) , (7) 

and h^jy denotes the regularized metric perturbations. This equation can be also reinterpreted 
as a geodesic equation for a point particle in the perturbed geometry g aj3 + h^p [26]. In order to 
evolve the SCO's trajectory subject to the self-force in a self-consistent way, we use the method 
proposed in [27], which is a relativistic extension of the osculating orbit method. The main 
idea is to model the SCO's worldline as a sequence of geodesies tangent to the worldline at any 
given point. The transition is smooth since EMRI trajectories remain close to a geodesic for a 
significant amount of time. In principle, to carry out this transition we need to take into account 
all the parameters that characterize the instantaneous geodesic (the principal orbital elements) 
and the position of the SCO within that geodesic (the positional orbital elements). In the initial 
version of the Chimera scheme [19] we only consider the dissipative effects of the self-force, 
i.e. those that only affect to principal orbital elements. Then, the evolution of (E, L Z ,C/Q) is 
given by 

^ = _At) a a dLz = M) a m 
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where the SCO self-acceleration o" F is related to the self- force by a" F = 1 F" F , and where Cc? 
and C,^ are the timelike and axial Killing vectors respectively, and is the Killing tensor. 

2.3. Multipolar Post- Minkowskian Self- Acceleration 

Since the self-force is not known for generic Kerr orbits, we approximate it here via a multipolar, 
post-Minkowskian expansion (see e.g. [20, 28, 29J). We begin by recasting the post-Minkowskian 



metric perturbation in the far field as the Kerr metric of the MBH g^ H in harmonic coordinates 
(we shall discuss the issue of coordinates in next section), plus perturbations induced by the 
SCO. The latter can be decomposed into conservative (time-symmetric) and dissipative (time- 
asymmetric) perturbations h^. Neglecting the former, we then have 

C = g^ H + C (io) 

The metric perturbations hf^: are the quantities we use to describe the radiation-reaction effects, 
and hence they are our approximation to the regularized metric perturbations h^p of Eqs. ^ 
and ([7]). The rate of change of the principal orbital elements in then described by Eq. ([§]), with 
the force of Eq. and the regularized metric perturbation of Eq. (10). 



The regular metric perturbation is a function of certain time-asymmetric scalar, V^ R , and 
vector, V£ K , radiation-reaction potentials, which to lowest order can be written as 

hft = 2 V m , = -4 K R , hff = 2 V m % , (11) 



1 For a point mass, the first-order metric perturbation is singular at the particle location, and a regularization 
prescription [25] must be employed to separate out the piece responsible of the gravitational backreaction. 



and the expressions for these potentials themselves are [2U1 ED] 



V RR (t«, * B ) = - \^ M iP (*h) + i^4 ife Mg(t H ) - i^M^^) + 0(x§AfW ) , (12) 
V^h, * B ) = 2l^< yfe >M|)(t H ) - ± eijk x%S$(t K ) + O(xTm^) , (13) 

where (xg) = (t H , x^) are spacetime harmonic coordinates, e ijfc is the antisymmetric Levi-Civita 
symbol, and 

x < ijk> ee xf - -xl8^x k J . (14) 
5 

is the symmetric trace-free (STF) projection of the multi-index quantity x Jjffc = x l x^x k . The 

quantities M^f , Afg, and are the nth-time-derivative of the STF mass quadrupole, mass 
octopole and current quadrupole multipole moments of the source. The first term in V RR 



[Eq. (12)] corresponds to the well-known Burke-Thorne radiation reaction potential |31j . 

The Chimera local self-force is completely specified once we prescribe how the multipole 
moments depend on the orbital trajectories. This can be achieved by asymptotically matching 
a PN and a post-Minkowskian solution in some buffer zone where both approximations are 
valid. The current Chimera implementation employs the leading-order expressions for these 
multipoles. As such, we are not consistently keeping a given PN order, but instead using 
leading-order expressions for the multipole moments, while keeping several such moments in the 
expansions. To lowest order, the mass moments are given by 



M tj = nm Z<ij> , M ljk = n5m z <ijk> (15) 

and the current quadrupole moment is 

Sij = r) 5m e kl<i Zj> k z l , (16) 

where r\ = m+M 9 /(m+ + M,) 2 is the symmetric mass ratio and 5m = M m — m^ is the mass 
difference. To higher order, these moments become more complicated as there are now non- 
linear contributions from the non-reactive potentials (i.e. non-linear contributions from the 
background) as well as tail and memory contributions. 



2-4- From Boyer-Lindquist to Harmonic Coordinates 

From the previous discussion it is clear that the MBH metric must be written in harmonic 
coordinates in order for all terms to be in the same coordinate system and to be consistent 
with other ingredients of the Chimera scheme, in particular the waveform generation procedure 
that we describe in next section. Nevertheless, there are parts of the scheme that are easier to 
implement in BL coordinates, such as the integration of the geodesic equations. Therefore, it is 
crucial to find a map from BL to harmonic coordinates. 

Harmonic coordinates refer to coordinate systems, {x^}, that satisfy Dig = 0, where □ is 
the D'Alembertian operator: □ = g Q ^V Q V^. In based on work by Ding (see Ref. [32] for 
references on this), a transformation to harmonic coordinates was found. Following [32], we can 
map the Kerr metric from BL coordinates (t,r,6,cp) to harmonic ones (£ H , x H , y H , z H ), and the 
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S.5. Waveform Generation 

Given the orbital evolution in harmonic coordinates one can construct the gravitational 
waveforms through multipolar, post-Minkowskian expressions in terms of a sum over multipole 
moments (see e.g. |33j for a review). This prescription is fully specified once expressions for the 
radiative moments are given in terms of derivatives of the orbital trajectory. Such identification, 
however, is difficult, as these moments are defined in the far-zone and have no knowledge of 
the source multipole moments, defined in the near zone. In this initial version of the Chimera 
scheme, we only consider the leading-order contributions to the multipoles moments, i.e. we 
identify the source and radiative moments. Then, in a transverse-traceless gauge and keeping 
multipole moments that include both the mass hexadecapole and the current octopole (thus 
keeping contributions one order higher than traditional kludge waveforms), we have: 
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The expressions for these multipole moments has been given in Eqs. (15) and (16), except for 
the mass hexadecapole and current octopole multipoles, which are given respectively by 
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The plus and cross-polarized projections can then be constructed via 



(20) 
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where e+ x are the plus- and cross-polarization tensors. Finally, the observables, i.e. the GW 
response function, are given by a projection of the plus- and cross-polarized waveform with the 
beam pattern functions of the detector. For a detector like LISA, there are two such functions, 
F+j/n and F x I /jj (see, e.g. [SI HI) and the response is thus 
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where D L is the luminosity distance from the source to the observer and the prefactor of v3/2 
is due to the triangular arrangement of the LISA spacecraft constellation. 



3. Numerical Implementation and some Numerical Results 

One of the main ingredients of the scheme is the integration of the geodesic equations ([2|-([5]). 
The resulting equations for the radial and polar BL coordinates, Eqs. ^ and Q respectively, 
have turning points (at pericenter and apocenter in the case of the radial coordinate, and at the 
location of the orbital inclination angle in the case of the polar coordinate). At these points we 
have either r = or 6 = and the ODE solvers may encounter convergence problems. To avoid 
these problems, we use angle variables associated with radial and polar BL coordinates via 



COS 9 = COS Oram COS \ ■ (23) 

1 + e cos ip 

For convenience, we parametrize the trajectory in terms of the BL coordinate time t, which is 
also a time harmonic coordinate, instead of the proper time r. Then, the direct result of the 
ODE integration is the sequence (t,tp(t),x(t),<f>{t)) for t £ [t ini , t enA ]. 

The main challenge in the numerical implementation of the Chimera scheme is the evaluation 
of the time derivatives of the different quantities involved. To understand this problem let us 
focus on the computation of our multipolar, post-Minkowskian self-force. The radiation reaction 



potentials V RR and V RR in (12) and (13) depend on up to the seventh time derivative of the mass 
quadrupole and octopole moments and up to the fifth time derivative of the current quadrupole 
moment. But since one also needs certain time derivatives of these potentials [see Eq. Q] one 
must compute up to the eighth time derivative of the mass quadrupole and octopole moments 
and the sixth time derivative of the current quadrupole moment. 

In principle, the computation of these derivatives can be done analytically by using repeatedly 
the equations of motion, Eqs. ([2])-([5]). The problem is that we need the time derivatives of 
the trajectory in harmonic coordinates, and to pass from the ODE angles (ip(t),x(t), 4>{t)) to 



harmonic coordinates we first need to use Eq. (23) to go from these angles to BL coordinates, 



and then Eqs. (17) to go from BL to harmonic coordinates. Therefore, we would need to obtain 
analytically higher time-derivatives of the ODE angles (up to eighth order), which involves using 
Christoffel symbols of the Kerr metric and several of their derivatives, and also to differentiate 



several times Eqs. (23) and (17). In practice, this makes the analytical computations unfeasible, 
even using modern computer algebra systems. 

Therefore, one needs to resort to numerical evaluations of these derivatives. First, notice that 
the trajectory, the velocity, and the acceleration can be computed with high accuracy directly 
from the integration of the ODEs in Eqs. Then, using analytical expressions, we can 

directly obtain the second time derivatives of the mass quadrupole and octopole moments and 
the first time derivative of the current quadrupole moment. Thus, we just need to compute 
numerically up to six additional time derivative of the mass moments (i.e. starting from their 
second time derivative) and up to five additional time derivative of the current quadrupole 
moment (i.e. starting from its first time derivative). 

Computing numerical derivatives is a subtle task, which we tried to implement in different 
ways. For instance, we experimented with many finite difference rules (involving from a few to 
more than 20 evaluation points), as well as other generic numerical differentiation techniques 
(such as numerical interpolation or Chebyshev differentiation). We found that, in order to 
obtain decent performance, a substantial amount of fine-tuning was necessary. To overcome 
this difficulty, we looked for a method better adapted to our problem. Notice that multipole 
moments are functionals of the trajectories, which are piecewise Kerr bounded geodesic, and 
bounded geodesies can be characterized by three fundamental frequencies (in the generic case). 
Then, an arbitrary functional of Kerr orbits, say f[tp, x, </>](*)) can be expanded in a multiple 
Fourier series of these frequencies as [23j 

m X ,m) = E 4, m ,ne- l( ^ +m ^ +n ^ }t , (24) 

k,m,n 



where (k, m, n) are integers running from — oo to +00 and fu mn are complex coefficients such 
that f_ k _ m _ n = fkmn- There are three special cases in which this expansion is simplified: 
(i) Circular equatorial orbits; (ii) Equatorial non-circular orbits; (iii) Circular non-equatorial 
orbits. In case (i), the Fourier series contains only a single frequency, the azimuthal one O In 
case (ii) and (iii), there are two independent frequencies, (Q r ,fL) and respectively. 

By modeling the required multipole moments as an expansion of the form of Eq. ( 24 ) , using 
a standard least-square fitting algorithm, time derivatives are then computed via 



fW[i>,X,<f>](t) = £ /i£»,ne" i(MV+mn ' +nn * ) *. (25) 
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where f k N m n = (— i) (Ml r + mf! fi + nQ ( j ) ) N f kmn - We tested this technique with different types 
of functions and we have found that it is very robust and provides very high accuracy even for 
the highest derivatives. For instance, the sixth time derivative is typically accurate to one part 
in 10 5 , which is more than enough for our purposes. 

In order to illustrate the strengths of the Chimera scheme, we present results for a generic, 
eccentric and inclined orbit. The evolution of all the constants of motion/orbital elements for 
many orbital periods is shown in Figure. [TJ While the large scale evolution looks very smooth, 
the shorter scale evolution in the subplots shows patterns with periods that match the orbital 
ones. These patterns are a consequence of using a local in time self-force and do not appear 
in evolutionary schemes based on flux averaging over a certain number of orbital periods. We 
can also appreciate from the global evolution that all quantities decay in time except for the 
inclination, which grows. However, if we look at the evolution over a few orbital periods, e.g. for 
the eccentricity, we can see that it can grow locally in time although the global tendency is to 
decay. The Chimera scheme leads to rich evolutionary patterns due to its local-in-time character, 
which also makes it a very valuable tool to investigate questions like transient resonances |21j . 




Figure 1. Evolution (for a total time of 10~ 2 yrs) of an eccentric and inclined inspiral of a system 
characterized by: M. = 1O 6 M , a/M. = 0.98, and q = 1CT 5 . The plots show the evolution of 
the following quantities: Energy E (top left), actually 10 A [E — E(t Q )], where E(t a ) = 0.9575513; 
angular momentum along the spin axis, L z (top centre); Carter constant C (top right); semi- 
latus rectum, p (bottom left), with p(t Q ) = 7; eccentricity, e (bottom centre), with e(t ) = 0.6; 
and inclination angle l (bottom right), actually 10 4 [t — t(t a )], where i(t Q ) = 57.39 deg. All plots 
contain subplots where the detailed evolution during a few orbital periods is shown. 



4. Conclusions and Discussion 

We have presented the main ingredients of the Chimera scheme, a framework initially designed 
for the modeling of EMRI dynamics and GW emission, but in principle also adaptable to 
intermediate-mass ratio systems. This scheme combines ingredients from the multipolar, post- 
Minkowskian formalism and black hole perturbation theory to evolve a non-geodesic wordline 
(with respect to the MBH geometry) and to construct waveforms. The orbits are built as 
a sequence of local geodesies, whose orbital elements evolve according to a local self-force 
that we approximate via a multipolar, post-Minkowskian expansion. The leading-order term 
of the self-force corresponds to the well-known Burke-Thorne radiation-reaction potential. A 
crucial ingredient in this construction is the mapping from BL coordinates, which we use to 
locally integrate the orbits, to harmonic coordinates, required both by the multipolar, post- 
Minkowskian self-force and by the GW multipolar expansion. Once we have trajectories in 
harmonic coordinates, it is straightforward to build waveforms. 

The Chimera scheme can be improved in different ways. One possibility is to use more 
accurate, beyond leading-order, expressions for the different multipole moments. Although these 
PN corrections are in general not known for generic spinning binaries, they should certainly 
improve the accuracy of the Chimera evolutions for certain orbits. Another way to improve the 
scheme would be to introduce conservative corrections to the background, for instance as it is 
done in the EOB formalism [351 1361 137]. 

A more detailed and exhaustive validation of the Chimera scheme would include, as a first 
step, a comparison of the evolution of the constants of motion with those obtained by solving 
the Teukolsky equation [38, 39]. One has to be careful in doing this comparison given that the 
latter employs averages over several cycles, while the Chimera fluxes are computed locally at 
the SCO's location. Once the fluxes have been validated, one should compare the waveforms 
themselves. An overlap study would determine the level of agreement between them. 

The Chimera scheme might also be useful to model IMRIs and even systems with moderate 
mass ratios. Recent comparisons between self-force, PN, and numerical relativity computations 
have concluded that by replacing the mass ratio q with the symmetric mass ratio r], the self- force 
predictions compare quite well with numerical relativity and PN predictions in the comparable- 
mass range |4U|, 141] . These results may allow for a simpler description of IMRIs, which 
otherwise would require either long, full numerical computations or higher-order perturbative 
computations, or a combination of both. In the Chimera scheme, q — > r] is actually mostly 
built in already, since the mass ratio information enters through the definition of the multipole 
moments. For these, we already use general binary expressions, instead of effective one-body 
ones with mass qM 9 . 

The ability to compute approximate local quantities at the location of the SCO, in particular 
the multipolar, post-Minkowskian self-force, makes the Chimera scheme an interesting tool to 
study certain local behavior believed to exist in EMRIs. For example, Flanagan and Hinderer |21| 
have recently reported that certain rapid changes in orbital elements can arise for generic EMRIs 
when the orbital frequencies become commensurate. During these rapid changes, it has been 
suggested that EMRI waveforms may present glitches. Reliable quantitative predictions for these 
glitches are needed to assess the importance of this effect for gravitational wave astronomy. In 
this sense, we could compare the effect of these glitches in waveforms as computed from an 
averaged-scheme (such as the Teukolsky one) and a local one (such as the Chimera scheme). 
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